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■ , Abstract 

^ ■ We study /3-ensembles with Bn, Cn, and Dn eigenvalue measure and their relation 

with refined topological strings. Our results generalize the familiar connections between 
local topological strings and matrix models leading to An measure, and illustrate that 
all those classical eigenvalue ensembles, and their topological string counterparts, are 
related one to another via various deformations and specializations, quantum shifts 
and discrete quotients. We review the solution of the Gaussian models via Macdonald 
identities, and interpret them as conifold theories. The interpolation between the 
various models is plainly apparent in this case. For general polynomial potential, we 
calculate the partition function in the multi-cut phase in a perturbative fashion, beyond 
tree-level in the large- iV limit. The relation to refined topological string orientifolds on 
the corresponding local geometry is discussed along the way. 
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1 Introduction 

Eigenvalue ensembles with measure to a power of /3, El widely known just as /3- 
ensembles, and their relation to topological gauge and string theories have been studied 
extensively in recent years. The special instance /3 = 1 of generalized interest is the 
Dijkgraaf-Vafa relation [1] between matrix models, supersymmetric gauge theory and 
the topological string. In more recent times, the focus has shifted to the more general 
situation with arbitrary /3, which relates the eigenvalue ensembles to fi-deformed gauge 
theories, refined topological string theory |2l[3] and the AGT conjecture Here, the 

^A]^ denotes the finite Coxcter group, and /3 a positive real number. We review the definitions in 
section [21 
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equivariant parameters of the f2- deformation, the ensemble parameter /3 and the 
string couphng gg are related via [2], 



One may note that neither the matrix model nor the topological string at present 
knows a microscopic interpretation for the deformation parameter /3. Rather, the 
mutual agreement of results calculated with different schemes, the consistency of the 
space-time interpretation via BPS state/instanton counting, as well as the relation 
with the fi-deformed gauge theory, especially in the Nekrasov-Shatashvili limit [5] 
give confidence that one should view all these models as integral part of a larger 
interconnected web of theories, thereby in fact defining various notions of quantum 
geometry, such as that of [3]. 

The prototypical example for much of this is the Gaussian model, with quadratic 
potential for the eigenvalues and corresponding, respectively, to a deformed conifold 
target space, (refined) Chern-Simons theory [6], as well as the c = 1 non-critical string 
at radius i? = /3 [TJ |2]. This Gaussian model also serves as building block for more 
general backgrounds. 

The purpose of the present paper is to take this logic one step further, and to study 
the possible role played by eigenvalue ensembles with other finite group measures, 
specifically, -Bat, Cn, and D^. These models, which we will refer to as Macdonald 
ensembles, are rather natural, and easily defined, but have been less studied in the 
recent topological string/gauge theory literature. The ensembles at /3 = 1 appeared 
briefly in the context of the Dijkgraaf-Vafa relation to four- dimensional A/" = 1 gauge 
theories with SO/Sp gauge groups and adjoint matter, and the realization of these 
gauge theories as string theoretic orientifolds. Most closely related to the spirit of the 
present work are [8] and [9]. Due to the nature of the original DV conjecture, these 
studies were essentially confined to tree-level. One of the aims of this work is to study 
the Bm/Cn- and Dj^ eigenvalue ensembles with general (3 beyond tree-level in greater 
detailtl 

It is then natural to expect that the Bn, and Macdonald ensembles with general 
/3 7^ 1 are related to a refinement of topological string orientifolds, which was one of 

^Since the root systems of Bpf and Cjv differ only in the length of the roots, hence the Haar 
measures are identical up to an overall factor (see for instance ^U\), it will be sufficient for us to 
consider only the i?jv and Djv ensembles. 
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the original motivations for the present workjj In thinking about the various pictures, 
it is however important to remember that these eigenvalue ensembles at /3 = 1 are 
in general not identical to the usual SO and Sp matrix models. Rather, the latter 
models provide the microscopic realization of the An ensemble at /3 = 2 and /3 = 
1/2, respectively. They are dual to A/" = 1 SO / Sp gauge theory with matter in the 
symmetric/antisymmetric representation. In particular, the orientifold in the large- TV 
dual topological string side acts differently on the tree- level geometry [121 US] • 

On the other hand we have the realization, in the Gaussian model, of the /3- 
parameter as the radius of the circle for c = 1 non-critical string. There, the orbifold 
of the c = 1 CFT at the self-dual radius is indeed equivalent to the R = 2 circle 
theory. This connection suggests the existence of an entire new branch of topological 
string/matrix model dual pairs that connects up to the standard branch at /3 = 2. Our 
work suggests that this is where the Bn and Dn Macdonald ensembles fit in. 

Whereas the duality between Af = 2 U{N) gauge theory softly broken to A/" = 1, 
An eigenvalue ensemble with /3 = 1, and topological string theory on the (spectral 
curve) geometry at large has been discussed and checked exhaustively in many 
works, for general /3 much less is known. For Bn and Dn, even at /9 = 1, no higher 
genus check of the proposed duality between the eigenvalue ensemble and topological 
string orientifolds has been performed. The power of (3 plays a major role beyond tree- 
level, and hence one might hope to be able to learn something about refined topological 
string theory and orientifolds thereof along the way, which are expected to be related 
to these /3-ensembles in the large N limit. 

In fact, under which specific conditions the eigenvalue ensembles for /3 7^ 1 relate 
to refined topological string theory in the large limit has not been pointed out 
so far in the literature, even for the ordinary An measure. Some examples where 
such a relation holds in a non-trivial manner where reported in [2], [3] . Although the 
calculations of [3] were restricted to the cubic ensemble at 1-loop level, there is little 
doubt that the observed correspondence extends to all genus, at least for the cubic. 
On the other hand, the Chern-Simons matrix models studied in [13] appear to indicate 
that in general such a relation does not hold. Attempts to formulate a refined version 
of the remodeled B-model of [15] have also failed to our knowledge so far. 

Some of the problems with the general applicability of the An type /3-ensemble can 

■^Meanwhile, the refinement of topological string orientifolds has been studied, with a different 
perspective and motivation, by Aganagic and Schaeffer [llj . 
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be traced back all the way to tree-level, that is, to the dual spectral curve geometry of 
the ensemble. This is most clearly visible at hand of the remodeled B-model geome- 
tries of [16]: In general, the spectral curve of the eigenvalue ensemble differs from the 
usual B-model target space geometry of the dual topological string and has singular 
points. Singularities are a general indication that refinement, i.e., a deformation of 
the correspondence away from /3 = 1, will fail. Indeed, singularities in the B-model 
geometry could harbor blow-up modes, which spoil an invariant BPS state counting. 
The corresponding mirror statement is the well-known fact that in order to have a 
well-defined BPS state counting of left and right spin (and not just the index), the 
A-model/M-theory geometry should be rigid {i.e., have no complex structure deforma- 
tions). This leads us to a condition on an A^r type /3-ensemble to have a well-defined 
BPS state counting interpretation. Namely, the spectral curve has to be non-singular. 
In particular, this applies as well to ensembles with polynomial potentials, i.e., one has 
to fill all cuts to ensure that one has a well defined BPS index. Under this restriction, 
the duality of [2] has a chance to survive the /3- deformation in a quite general setting. 
Similar considerations apply to the and Dn cases, up to some technicality which 
we will explain in more detail in section 15. 1[ Confirmation for this expectation will be 
found at hand of A^r, Bj^ and -Djv /3-ensembles with quartic potential, which appear 
to be as well compatible with a (refined) topological string interpretation, as the free 
energies fulfill the 1-loop holomorphic anomaly equation. 

The outline is as follows. In section [2] we will give the definition of Macdonald 
ensembles with special emphasize on An, Bn and D^. This is followed by a detailed 
discussion of the large expansion of the Gaussian partition functions and implica- 
tions thereof for refined topological string orientifolds, in section [31 In section H] a 
recursion relation satisfied by Gaussian correlators is derived (generalizing [T71 [T5]). 
which constitute an essential ingredient for the explicit calculation of the multi-cut 
ensemble partition function, which section [5] is about. In subsection 15. ![ we will give a 
generalization of the framework of [19], [201 113 IH] to B^ / D^, and apply it in section 
15.21 to the model with quartic potential. The B-model verification of the tree-level and 
1-loop results of section 15.21 will be performed in section [SI We conclude in section [3 
In appendix [X] the explicit results for the Qs expansion of the free energy of the B^ 
and Dn /3-ensemble with quartic potential are attached. 
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2 Macdonald ensembles 



Let Q he a. finite group of isometries of generated by reflections in liyperplanes 
tlirougli tlie origin {i.e., a finite reflection or Coxeter group). Let tliere be h hyper- 
planes, each defined by a condition on A G of the form 

N 

where a„ G M^, a = 1,... , h. The group Q naturally acts on the algebra of polynomial 
functions on M^. The ^-invariant polynomials form an M-algebra generated by 
independent polynomials of degrees di, i = 1,. . . ,N. Normalizing the vectors {aa,i)i 
via X^ili ^a,i = 2, define the particular invariant polynomial 

h N 
0=1 i=l 

Macdonald conjectured the integral identity [21] 

where [d\] := YliLi dXi, (3 E C with Re/3 > 0. A proof of this identity has been given 
by Opdam [221 [23]. Zg{P) is also referred to as Macdonald integral. 

For Q = An^i the Macdonald integral specializes to Mehta's integral. For this 
work, in addition the cases Q = Bj^ and Q = Dj^ are of particular interest. These give 
rise to the integral identities. 



i=l i=l 

N-1 



ZdAP) - JJdX] A(A ) e ^ - ^(1 + ^3) r(l + /3) ' 

where A(A) denotes the usual Vandermonde determinant, A(A) := ]n!i<j('^j ~ -^j)- We 
are particularly interested in the large limit thereof, see section [31 

Viewing the above integrals as partition functions of Gaussian eigenvalue ensembles, 
it is natural to define general Macdonald ensembles by replacing the quadratic term 



6 



J2i ^1 with a general "single-trace" polynomial potential W{\i). For Q = Ajv_i the 
Macdonald ensemble is identical to the usual /3-ensemble. For Q = Bj^j and Z^jy? and 
with /3 = 1, these ensembles correspond to the ones considered in [8] in the context of 
the Dijkgraaf-Vafa relation (with the additional condition W{x) = W{—x)). 

It is convenient to parameterize the measure P(b^rf)(A) := Pg{X) with Q being A^, 
Bj^ or as 

N 



P(M)(A)=A+(A)^+'^A_(A)nA^ (2.3) 



where we defined 



N 



i<j 

In particular A_(A) = A(A), corresponds to the usual Vandermonde, and A_(A)A+(A) = 
A(A2). For {b,d) = (0,0) we get PaA^), (1,0) yields PBt^W and (0,1) results in 
Pdn{X). Hence the Macdonald ensembles Zg(/3_) with Q = An, Q = Bn ot Q = 
can be treated simultaneously via the ensemble l 

^{M)(/3) - j MA]|P(M)(A)r^e-^-^(^«) , (2.4) 

with P(6 (i)(A) as defined in (12. 3p . 

The expectation value for an operator insertion O is defined as usual as 

:= j MA]|P(M)(A)P'^ae"S-^(^'). 

Trivially, we have (l)(b^^) = Z(^b,d){f3). 

It is instructive to compare the formulas for Bn and D^- The only difference is the 
additional factor of YliLi ^1.^ for B^, and can be interpreted as follows. We know from 
the An /3-ensemble that the insertion of a brane at position x in the dual geometry 
corresponds to the insertion of some power of a determinant factor of Hilil^^Aj) times 
an overall classical piece of ^d{x) = e^^^^ [3]. Different powers of the insertion 
correspond to different types of branes [3]. Thus, the insertion of a brane at x plus a 
mirror brane at — x corresponds to an operator insertion of 

N 

^^{x)^f,{-x) = - ■ (2.5) 

1=1 



^The Qs dependence needed to match to topological strings will be brought in via W{x) and, if 
necessary, a rescaling of the eigenvalues. 
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This implies that the Bjsf ensemble can be understood as the Dj\f ensemble with inser- 
tion of an additional pair of branes at the origin (x = 0), i.e., 



Z(i,o)(/3) ~ (vi^;3(0)vi^^(0))(o,i) = Z(o,i)(/3) , (2.6) 



where we defined the partition function with h coincident (/?-) branes \E'|V^,Jx) in the 



background parameterized by (6, d) as 

((*^w)"> 



<.,w^- (2-7) 



3 Gaussian ensembles: Large N partition functions 

In this section, we review in some detail the large- iV expansions of the Gaussian par- 
tition functions (12. 2p . as well as the various ways that these enter into the topological 
string. 

3.1 Large expansions 

An 

In contrast to the and ensembles, the 't Hooft large- limit of Za^_^{P) 
has been studied extensively in the physics literature. The asymptotic expansion as 
— oo is related to the "Schwinger" integral, 

f dt e~^* 

logZA^_M- J 7 (eM-i)(e^2f-i) =-'^ogZA{g,,f3), (3.1) 
where ei, €2 are related to gg, (3 as in (II. ip . and we have 

/.:=V^^,iV. (3.2) 



We note the obvious symmetry of the partition function under ei O €2. As — t- 0, we 
have the well-known asymptotic expansion 



] t (e^^*-l)(e-*-l) ' ^"-"^ 

with certain polynomial expressions $^^(/3). [f| 



^Really, /3"$^-* is a polynomial in /3. 
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As is well-known, for n even, the ^^^\/3) specialize at /3 = 1 to give the virtual 
Euler characteristic of the moduli space Aig of genus g = ^ + 1 complex curves 



if--Hl)^x(M,)^^£^^, (3,4) 

where Bn are the Bernoulli numbers. For n odd, the $^■'(1) vanish. 

On the other hand, for n odd, the specialize at /3 = 2 to give the virtual 

Euler characteristic of the moduli space Ai'? of complex curves of genus g = n + 1 
with a fixed-point free anti-holomorphic involution [i.e., certain type of real curves) 
^2E\ [271 ES], up to a rescaling of gs- In string theory language, the quotients give 
unoriented Riemann surfaces with genus g = g/2, [g being even), no boundaries, and 
one crosscap. 

For n even, we have that $^^^^''(2) = 2~^x(A^g). Hence, the coefficients at /3 = 2 
show the typical structure of an orientifold 

2"M")(2) = ix(-M,)-x(-M,^). (3.6) 



For later reference, note that via making use of fl3.1l) and (13.61) . one can infer as well 
an integral representation of the generating function for the x(A^^); ^-e-? 

ZA{V2g,,2) 1 f dt e-^' ^ ( uo J dsy"' 



(3.7) 



A similar, and in fact related, "Schwinger" integral has appeared before in the context 
of the orientifold constant map contribution [291 EO] • 

The fact that the Aj^ /3-ensemble can be used to interpolate between the Euler 
characteristic of moduli spaces of complex and real curves was pointed out in [27] , and 
interpreted as a geometric parameterization. In particular, it was conjectured that 
the $^'*(/3) themselves should describe the Euler characteristic of some related moduli 
space. 

Although the appearance of the moduli of real curves is suggestive, a simple closed 
string theory interpretation is hampered by the fact that the expansion (13. 3 p for /3 7^ 1 
contains terms of both even and odd powers of gs- This originates from the fact that 
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(13.11) is not invariant under (ei, 62) — > (— ei, —€2), except when ei = —€2- As is by now 
well-appreciated, the additional (quantum) shift 

/.^/. + ^^ (3.8) 

restores that symmetry. We have the asymptotic expansion 

with \l/^''(/3) = for n odd. This shifted partition function is also identical to the 
partition function of the c = 1 string at radius R cc P, originally found in [21]. From 
the above formulas it is clear that $^■'(1) = \1/|^^(1). 

Turning to the topological string, it was discovered long time ago in [7], that the 
integral (13. ip at /3 = 1, ^.e., the c = 1 string at the self-dual radius, governs the leading 
behavior of the B-model topological string in the limit in which the target space devel- 
ops a conifold singularity, as the complex structure parameter yU — )■ 0. As explained for 
instance in |32], the coefficients $^■'(1) = \1/^''(1) therefore provide universal bound- 
ary condition for solving the topological string via holomorphic anomaly equation [33] . 
As shown in [311 EH], see also [36], the one-parameter deformation '^^2\f^) provides 
the analogous boundary conditions for solving the refined topological string in the B- 
model via the same holomorphic anomaly equation. (Alternatively, one may use the 
extended holomorphic anomaly equation of [37] with boundary conditions provided 
by the $^^(/3) to solve for the refined topological string amplitudes after undoing the 
quantum shift.) This observation confirms the identification of /3 as the radius R of 
c = 1 string [2] . 

A seemingly unrelated observation is the fact that the coefficients \l/^^(2) also have 
a topological string interpretation, in the context of the real topological string [38] . 
Namely, writing [34] . 

2n/2^(n)(2) = + X[/kb) , (3.10) 

the \I/kb control the leading behavior of the topological string amplitude on a genus g 
Klein-bottle (an unoriented Riemann surface with genus g and even number of cross- 
caps) around a conifold point in moduli space. This relation begs for a topological 
interpretation of the \&kb similar to that of the \1/^^(1) in terms of the moduli space 
of genus g complex curves (with n = 2g — 2). We note that whatever this interpre- 
tation is, \1/^^(1) is not the virtual Euler characteristic of moduli of complex curves 
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with fixed-point free anti-liolomorphic involution of odd genus g = 2g — 1 (which are 
the covers of these higher genus Klein bottles) studied in [27], which vanishes, but 
should be closely related to it. It is also interesting to note that the quantum shift 
(13. 8 p transforms Klein-bottle contributions into cross-cap contributions, as is apparent 
via comparing (13. 6p and (I3.10p . 

Bn 

It follows from (12. 2 p that Zbj^{P) can be expressed in terms of Zaj^_j^{/3) as 

r(l + 2/3) 



logZsA^) = log ZA,_,m + iVlog 



r(i + /3) ■ 



Since in our 't Hooft limit, we neglect the most singular terms, of positive power of A^, 
we simply write, 

Zs(^7„/3)~ZA(^?s,2/3). (3.11) 

As a result, the /3 -H- 1//3 symmetry of Za{P) translates to a /3 o 1/(4/3) symmetry 
of the Zb{/3) partition function. Similarly as in the Am case, we denote the expansion 
coefficients of the Qs expansion of the corresponding (shifted) free energy as \l/^''(/3). 
Obviously, 

4")(/3) = vl/W(2/3). 

For later reference, let us explicitly state the "1-loop" coefficient which the (shifted) 
ZBigs,P) implies, i.e., 

^2^/3) = ^ Q + 4/3) . (3.12) 

Dn 

For the Dn Macdonald ensemble, we sort terms such that 

logZ„„(,?/2) ^ logZ.„_.(« - ,og jni±|| H- A^log j|l±|L . (3.13) 

Thus, besides the An term, we have one additional dependent term which we don't 
neglect in the large- A^ limit. Invoking the integral representation of the digamma 
function 7(x) = £\og(T{x)), 



7(x) = / dti — 
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one can infer that the essential part of the new contribution reads 



t 1 - e-* 

Redefining t — t-^ and taking — oo while keeping y/^Ngs =: fj, fixed yields 

dt e-^* - e^/^'/^ 



t 1 _ e v^* 

Reverting to our usual notation, we obtain (up to non-universal terms) 

logZ,„(^/2) ~ /f + / f =: logZ,fc,3/2) . 

(3.14) 

We recognize the second term as being essentially the generating function (13. 7p . Hence, 
Zoids, P) is entirely given by a combination of Za, i-e., 

log ZDigs, /3) = log ZA{gs, 2/3) - 2 r(^(7.) . (3.15) 

(where T is defined in (13. 7p .) It is also instructive to express the "1-loop" coefficient 
$2^(/3) contained in flXTil) in terms of $i°H/3)- Using the relations fl3A5ll and flO) . 
we infer 

$g)(/3) = $i°)(2/3) + 2 ($5^(2) - ^$?^(1)) . (3.16) 

We observe that this result matches the structure of the orbifold branch partition 
function for the c = 1 CFT on the torus derived in [32] • We take this as a hint that 
(I3.15P is related to the partition function of the orbifold branch of the c = 1 non-critical 
string. More precisely, such a relation should hold after an appropriate quantum shift. 
One way to identify the appropriate Z^Ar-analog of (13. 8 p is to impose a symmetry under 
gs —7- —gs- The symmetry can be motivated as follows. We know that for integer values 
of /3 the partition function of the c = 1 string on the circle branch can be matched to 
the partition function of the topological string expanded near a, Ap_i type singularity 
(under appropriate choice of deformation parameters) [ID]. Since the chiral ground ring 
manifold of the c=l string on the orbifold branch corresponds to a Kleinian singularity 
of -D-type [H], we expect that similarly the orbifold branch partition function can be 
matched to the topological string expanded near a D-type singularity, implying the 
symmetry under gs —gs- 

Indeed, after 

(3.17) 
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the free energy (13.141) becomes 

r ri+ e"''* cosh { ^^i±^2)^^ . . n 

/f 2sl.hM/2)s.:h(4 -i:-^i;'W<f) . (3.18) 

and clearly possesses an even power only expansion in gf^. § We denote the expansion 
coefficients by \l/^^(/3). For the special value /3 = 1, we find in addition 

^(^)(i) = = m^2\2) . (3.19) 

This identity is precisely the one expected from the matching of circle and orbifold 
branches of c = 1 string at i? = 1 and R = 2, respectively. It is important to note 
however that generally, Z^^ and Z^^ are not related by a simple shift of A^. 
For later reference, we explicitly state the "1-loop" coefficient 

3.2 Implications for toric Calabi-Yau backgrounds 

The free energy J^a{Q] P) of the refined topological string on the resolved conifold 
geometry, i.e., 0{—l) © 0{ — l) — )■ P^, is related to the refined deformed conifold free 
energy given by the integral (13. 9 p by identifying the Kahler parameter as Q = and 
simply replacing the integral by a sum 

/oo 
dt^^. (3.21) 
d=i 

This replacement originates from the sum over states of DO-brane charge A; G Z and 
mass ~ yU + 27rik, or, in M-theory language, from the extra state degeneracy due to 
momenta around the M-theory circle. It is natural to assume that a similar "quanti- 
zation" as in (I3.2ip can be applied to the other Macdonald integrals as well. For B^, 
we infer from (13. lip 

J^BiQ;P) = J'A{Q;2f3), (3.22) 

i.e., the resolved conifold free energy of type agrees with that of type A^. 

®It is interesting to note that this generating function looks very similar to the generating function 
for the massless hypermultiplet contribution occuring in SU(2) gauge theory on fi-deformed Ai ALE 
space [42] . 
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The Djsi ensemble is more interesting. Applying fl3.2ip to the (shifted) free 
energy fIS.lSp . we obtain 



d=l 



with the usual definitions q := e^^" and t := , as prediction for the refined free 
energy oi Dj^ type of the resolved conifold geometry. In particular, we have an even 
power only expansion in Qs and the relation J^oiQ] 1) = J^a{Q] 2) holds. 

It is instructive to compare this result with the expectations based on a topological 
string orientifold interpretation. A convenient reference is the recent proposal [TT] for 
the orientifolded and refined resolved conifold free energy. This proposal, which is 
obtained from an S0{2N) refined Cher n- Simons/geometric transition point of view, 
reads, 

d=l 

(we have here exchanged q ^ t (corresponding to /3 4-7- 1 //?).) As observed in [TT], the 
specialization of f l3.24p to /3 = 1 equals the free energy of an orientifold of the resolved 
conifold (acting either in fixed-point free |13], or in a real [HI [29] fashion). The second 
term in 03.241) can be understood as originating from the second term in (13. 7p of the 
(unshifted) free energy, summing only over even DO brane charge (up to a shift). It 
may also be seen as a brane placed at —1/2 log Q in the A-model geometry. Since the 
brane is localized in two space-time dimensions, it is exposed only to a single parameter 
of the f2-deformation (after a suitable redefinition of parameters). 

The structure of the refined orientifold free energy fl3.24p is consistent with the 
results of [31] , where it was found that the free energy of the fixed-point free orientifold 
of C(— 2) ©(9(— 2) — )■ X equals the refined free energy at /3 = 2 (for this orientifold 
one has no open string sector). 

We draw attention to the fact that, at /3 = 1, the shift of logQ in the first term 
of fl3.24p cancels the part of the sum in the second term coming from even d. That 
summation would be unusual for the Qs-odd sector of an orientifold. The cancellation 
is possible because the open string contribution is essentially a closed string period. 
On the other hand, this kind of comparison challenges the extrapolation of the pro- 
posed orientifold structure fl3.24p to more general toric Calabi-Yau geometries, since a 
cancellation between a shift and an open string contribution is not possible in general. 



14 



Instead, we propose to view J^oiQ] P) of (13.231) as the refined free energy of type D]^ 
of the resolved conifold, independent of an orientifold interpretation. Although for 
specific situations (such as the conifold or Dijkgraaf-Vafa type geometries, see section 
[5]), when the orientifold contribution is a closed period, the (unshifted) free energy 
can be matched and interpreted at /3 = 1 as an orientifold, this relation is not expected 
to persist in general. 

Specifically, we propose to identify the theory oi type with the orbifold branch 
in the c = 1 moduli space, extending the identification of An with the circle branch of 
c = 1 [2]. The relation P = R is the same on both branches. For a toric Calabi-Yau 
manifold, one may then use (I3.18P as boundary condition on the holomorphic anomaly 
equation at the conifold point in moduli space in order to obtain predictions for the 
theory. We will not pursue this quite interesting toric direction further in this work, 
but rather stick to the Macdonald ensemble setting, where the correspondence with 
orientifolds holds provided we work with even potentials. This will give further evidence 
via explicit calculations that the deformation (12. 4p away from /3 = 1 is consistent. 

4 Gaussian ensembles: Correlators 

In order to evaluate perturbatively /3-ensembles with multi-cut support we will need 
to evaluate normalized Gaussian correlators defined as 

c££ ,M-- '"y"''''-*' =^ /l''A||p.M,(A)pvia.^;n5,.. (4,1) 

with Sk := Xlili normalized via the Gaussian partition function Z(^h,d} given in 

iOl- Clearly, 

cK.,o(/3) = N X N X ■■■ X N . 

Correlators with non-vanishing ki can be solved for recursively invoking the Ward- 
identities resulting from invariance under eigenvalue reparameterizations. For example, 
more recently this approach has been followed in the An case with general /3 in |18j . 
The generalization of this approach to the Bn and Dn case we are interested in is 

^It is conceivable that this resuh (instead of p.24p ) can be obtained also from the Chern-Simons 
point of view by appropriately incorporating the quantum shift p.l7[) in the large- TV limit. 
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straight-forward. The Ward identities read 

g I [d\]dx, (^\- (^A+i\r'^A4\)f[\^^ e-^S-^?5fc,...5,„ I =0. (4.2) 

Acting with the derivative on each factor and expressing the resulting terms through 
the correlators fl4.1l) yields a recursive equation for them. For brevity, we will here 
explicitly state only the new contributions of A_|_(A)''"'"'^ and Hill -^i which do not 
occur in the An case. The remaining An contributions can be deduced similarly. 
Using, 

k=l \i=l J \i=l 

the new contribution which only occurs for Bn can be inferred to be simply given by 
The derivation of the contribution of A_,„(A)^+^ goes as follows. One rewrites 

N N 

J2Kdx, (A+(A)2/'(^+'^)) =2/3(6 + rf)A+(A)^^(''+'^)5^A^9,,logA+(A), 

k=l k=l 

and deduces the identity 
2/3ib + d)Y,Kdx,'^ogA^{X) = {b + d)/3Y,j-^ = ib + d)PY. A. - (-A ) " 

fc=l i^j * i^j * ^ 

Note that the second step in the above equality is only valid for n odd. Fortunately, 
knowing only how to deal with the n odd case is sufficient for our purposes. Finally, 
making use of the identity 

\n _ {_\ \n '^-i 

i^j * ^ k=0 ij^j k=0 

one arrives at the contribution 

n.-l 
k=0 
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Combining all terms, taking also the usual contribution into account (see for in- 
stance [HI US]), the Ward identities (14. 2p translate to the recursive relation 

n— 1 m 
k=0 j=l 

(4.3) 

Note that the sole difference between the S^v and Z^at case is a switch of sign of one of 
the terms entering the recursive relation. 

As pointed out already above, in case oi Bjq and D^^ the recursive relation is only 
valid for n odd. However, the relation closes if furthermore all ki are even. Hence, 
it can be used to determine all correlators C^^'^^ fcm "^i^h ki even. A few examples 
follow. 

Cp^^ = N{1 + 2P{N - 1)) , Cp^^ = N{1 + 2f3N) , 

d'f = Cr\2 + Cr') , cSf = Cp'\2 + C^°)) , (4.4) 
^(0,1) ^ ^(0,1)^3 ^ 4^^^ _ ^ ^(1,0) ^ '°^(3 + /3(4iV - 2)) . 

One should note that while in the case with /3 = 1 a generating function for the 
1-point correlators (7^"'°^ is known [25] and there are also closed expressions for (3 = 1/2 
and (3 = 2 [16], no such closed formula has been found for general /3, nor for the Bn and 
IPtv cases, so far. Nevertheless, we can make at least one general observation regarding 
the structure of the Cn'''^ for n even. Namely, the coefficients of the highest powers in 
appear to be always given in terms of the Catalan- numbers C„ := ; ^-c, 

c£'°H/3) = C„/3"Ar«+i ^ _ _ _ ^ c£°)(/3) = C„/3"A^"+^ + . . . , 
C(0'i)(/3) = 2"C„/3"A^"+^ + ... . 

5 Multi-cut potentials: Perturbative calculation 



(4.5) 



In this section we shall evaluate the eigenvalue ensembles (12.41) in a perturbative fashion 
via a saddle-point approximation, making use of the fact that for (7^ — )■ the eigenvalues 
Aj localize to the critical points of the potential W{x). A detailed exposition of the 
perturbative calculation of the /3-deformed A^ partition function has been given in 
|18j . as a generalization of the earlier /3 = 1 works [191 |20l [17]. We will not repeat that 
discussion here, but focus on the new features that appear for Bjsf and ensembles. 
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We discuss only potentials with the symmetry W{x) = W{—x). In particular, W{x) is 
a polynomial of even degree. 

5.1 Saddle-point approximation 

Since the degree d of W{x) is even, we have an odd number c = d — 1 of critical 
points. In particular, due to the Z2 symmetry of the potential there is one critical 
point at Im(x) = Re(x) = 0. We denote the set of critical points as /i*^'^-* with k E 
{ — 0, ^}, hence /i^"'^^ = The saddle-point approximation requires 

us to distribute the eigenvalues between the c critical points. Let us denote the 
number of eigenvalues located around yU^'^^ as A^^. In contrast to the Ajq case, we 
have to impose some additional constraints onto the eigenvalue distribution. This will 
allow the -Bat and Dn ensembles to be dual to, both, Af = 2 SO/Sp gauge theories 
with adjoint broken to A/" = 1 by a tree-level potential of the form W{x) and to 
topological string orientifolds. From an orientifold point of view, it is more convenient 
to work in the quotient space perspective. In particular, for the eigenvalue ensembles 
this allows to avoid to deal with the interactions between "mirror" eigenvalues under 
the Z2 identification of cuts which the duality to orientifolds requires. Following [H], 
we implement the quotient into the eigenvalue ensemble by localizing the eigenvalues 
around the quotient set of critical points, i.e., we take the eigenvalues to be localized 
around (c + l)/2 of the critical points such that 



The Aj\f condition for consistency of the /3- deformation of filling all cuts stated in the 
introduction then changes in the Bn/Dn case to filling only the quotient set of cuts. 

(k) 

The partition function can then be evaluated by considering small fluctuations i/n , 
with /c G {0, . . . , around the critical points, i.e., we set 



N = No + ... + Nc^ . 



(Ai,A: 



'2, • • • 




) 
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Under this decomposition we have that 



A_ (A) ^ n n - yn n n n (/^^'"^ - /^^"^ + yi^^' - yf 

k=0 i<j o<m<n<^ «=1 J = l 

c-l 



2 JVfc 



^4^) ^ n n (2/^^^^^ + ^J'' + yf) n n n (^^'"^ + ^^"^ + + 

(5.1) 

Hence, 

with interaction term 

oo I 



(5.3) 

The potential decomposes as 



i=l k=0 71=1 



(5.4) 



Note that the Bn, respectively ensemble decomposes in the saddle-point approxi- 
mation into (c — 1)/2 An and a single Bn, respectively Dn eigenvalue ensembles, which 
are coupled via the interaction term X(?/), as expected. Hence, after an appropriate Qs 
dependent redefinition of Sl"^\ and expansion in Qs to bring down powers of ^J'"^ the 
partition function of the eigenvalue ensemble (12. 4p reduces in the saddle-point approx- 
imation to a sum over Gaussian correlators which are determinable via the results of 
section m Due to the normalization of the Gaussian correlators, the resulting partition 
function has to be supplemented for each cut by a factor of Zg{P) (defined in (12. ip ) 
with Q either B^ or for the fixed cut and else An- 



19 



5.2 Examples: Bjy and Dj\f quart ic 

Let us now consider an example in more detail. The simplest non-trivial Z2 symmetric 
example is given by the quartic potential 



W{X) = -\' - . (5.5) 



Clearly, W{—X) = W{\) and the set of critical points jj^''^ consists of 

hence possesses a three-cut structure. The partition functions Z(^b^d) for A^, Bn and 
measure can be obtained in a perturbative fashion as outlined in the previous section. 
While one has to fill in the An case all three cuts as discussed in the introduction, for 
Bjy and Djy one has to fill only two of the cuts in order to incorporate the Z2 quotient 
as mentioned in section 15.11 

For the A^ case, let us just quote the relevant observation without giving any 
further details. Namely, we observe that if we fill all three-cuts, the disk sector {g~^) 
of the corresponding free energy is a combination of closed periods, i.e., 



(1/2) 




where as usual Si := NiQ^. This indicates that one has to perform the additional 
quantum shifts 

S^■.= (^^-^1-^))^- (5.6) 

in the large limit. Indeed, after the shifts, one obtains an expansion into only 
even powers of Qs of the free energy J-"^, as is necessary for a well-defined BPS index 
interpretation of the corresponding partition function. 

Let us now consider the B^ and Dn cases. For the quartic potential (15. 5p . the 
measure specializes to 



N 

b+d ■ 



i=l 
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the interaction term (15.31) reads 



Xiy) = - 2/3 f: ^ ^ Q (1 + (6 + d)i-ir) 

1=1 r=0 ^ ^ 

Z=l r=0 ^ ^ (=1 

and the potential contribution fl5.4p is given by 

oo _ 

^(^) = - E - {(d'^Wmsl^^ + {d-W){6)Sl^^) . 

Tl. 

n=l 

After performing the rescalings 

o{0) I __9^\ oil) ^ 3s \ ^(1) 



const. , 

(5.7) 



the partition functions can be expanded in and reduce to a sum over Gaussian 
correlators, which one can efficiently calculate following section HI For the reader's 
convenience, the explicit free energies to some lower order in gs are given in appendix 
Rl Defining 

So = 2Nogs, Si = N,gs, (5.8) 

we obtain 

^b\^Oi Si) = J^^^\So, Si) = -T^\So, Si, S^i = Si) , (5.9) 

which is the expected tree-level result from a topological string orientifold perspective. 
The first order open string corrections (g^^) read 



N dSo 2 V /3; dSi ^g^^Q^ 



(1/2) _ 1 d^S^ ^ 1 / i\ 



2/3 dSo 2 V dSi 



and are combinations of closed string periods. Note that for /3 = 1 we have that J^^^'^^ = 
_ _iQ^^jr^^s^^s^^g_^ = g^^^ con&ming the earlier resuhs of [H SH HHl [9] . 
In the dual topological string orientifold, the sign difference translates into the two 
possible choices of charge of the orientifold fixed-plane. Similar as for the An case, the 
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order ^ given in flS.lOp suggests to perform the additional quantum shifts as in (15. 6p 
such that 

T-(i/2) _ _ 1 dTf\So, Si, S^i = Si) ^(1/2) _ 1 d^^\So, Si, S-i = Si) /c- 1 1 n 
- 2 dSo ' "2 dSo ' ^ ^ 

and the relation J^^^"^^ = —T^j^'^^ continues to hold under the /3-deformation. In 
particular, the open string contribution at order is independent of /3. However, for 
higher orders in Qs, one has that generally 

and equality (up to overall sign) only for /3 = 1. Since the open string contribution is 
trivial {i.e., it is a closed period), it is more convenient to shift away the complete open 
string contribution such that J^^^'^^ = T^^^'^^ = 0, for g odd. Hence, in this specific 
gauge the free energies possess an expansion into even powers of gs only. If necessary, 
the open string contribution can be easily reinstated via performing a reverse shift. 
The main advantage of this shift is that it allows us to utilize the usual holomorphic 
anomaly of [33] instead of the extended holomorphic of [3Z] to reproduce the -Ba? and 
Dn partition functions in the B- model. On a technical level, the former is easier to deal 
with. However, we like to stress that the latter is more general, since it is expected 
to capture the partition function independent of any shift of parameters, similar as 
observed at hand of gauge theory in [311 ES] . 

Finally, let us comment on the Nekrasov-Shatashvili limit of the free energies (in 
the gauge with an even power gs expansion). In our parameterization the limit of [5] 
corresponds to 

>Vi^^= lim /3f ji^V/3) . (5.12) 
From our explicit computations we observe that 

i>v(f) = >vt^ = Wl^^ . (5.13) 

A similar non-uniqueness property of the limit has been already observed at hand of 
gauge theory on ALE space in [12], and is in fact as expected. This is because, since 
gs ^ (in order to keep h := fixed, cf., (11. ip ). the Nekrasov-Shatashvili limit, and 
hence W^^\ is intrinsically of tree-level nature. More specifically, the limit corresponds 
to a (semi-classical limit of a) quantization of the spectral curve of the respective 
eigenvalue ensemble, following [3]. Since the spectral curves of the A]\r, B]^ and D]^ 
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ensembles (under a proper Z2 identification of tlie A^r spectral curve) are identical (c/., 
(15.91) ). so should be the quantization thereof. The relation (I5.13P shows that this is 
indeed the case. 



6 B-model verification of and quartic 

6.1 Tree-level geometry 

The dual tree-level geometry of the quartic eigenvalue ensemble with S^v and D^^ 
measure with /3 = 1 has been discussed already in the literature to some extent (c/., 
[IHIEIIS])- Since the power of (3 is only relevant at one-loop and beyond, we essentially 
can borrow the known tree-level results. 

The periods of the dual geometry of the eigenvalue ensemble with potential W{x) 
of the form (15.51) are given in terms of the periods of the hyperelliptic curve 

y = Mix)^M^, (6.1) 

with 

M{x) = g , a{x) = \ {W\xf + /(x)) , (6.2) 

g 

and where f{x) is a degree two polynomial. In particular, the moment function M{x) 
is a constant because we fill all cuts. For the quartic, the curve (16. ip is of genus two. 
The effective one-form of the dual geometry reads 

cu = ydx , (6.3) 

which we express in terms of the six branch points Xi of the curve as 



w = 5-, 



6 



\ i=i 



The cuts are chosen to be [xi, ^2] , [^3, X4] and [^5, Xq] on the real axis. Imposing the Z2 
symmetry under x — ?• —x of the x-plane (this requires that f{x) is even, i.e., f{x) = 
62^^ + bo with bi parameterizing the complex structure) leads to the identification 
Xi -H- —Xq, X2 ^ —X5 and 0:3 ^ — X4, yielding the one-form 

u = g\J (x^ — xf)(x^ — xl){x'^ — x|) dx . 

The Z2 symmetric x-plane of the geometry is illustrated in figure [H 
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V 



Figure 1: The symmetric x-plane of the quartic with cuts, period contours and reflection 
symmetry indicated (covering space perspective). 

One should note that if one adjusts such that 62 = 0, the curve (16. ip takes the 
form of the Seiberg-Witten curve of four dimensional M = 2 SU (3) gauge theory, with 
choice of Coulomb parameters Oi = —02 and 03 = (the zeroes of W'{x) correspond 
to the a,j), which can also be matched to the curve of Sp{2). It would be interesting to 
extract the gauge theory gauge coupling and 1-loop gravitational correction following 
[50| [20] and see if one can match to a f2-deformed gauge theory, as is the case for the 
cubic and f2-deformed SU{2) [3]. However, since 03 = 0, we are not on the Coloumb 
branch, and things are expected to be a bit more tricky. Therefore we will not follow 
this rather interesting direction further in this work. 

Comparison with the one-form ( 16. 3 p expressed via ( 16. 2 p gives the relation 

S' = l{xl + xl + xl) . (6.4) 
In order to explicitly calculate the period integrals, it is useful to change variables to 

Zo = X3, 21 = ^(X2-Xi), X=^(xi + X2). (6.5) 

Hence, we set 

X3 = Zq, Xi=I - Zi, X2=I + Zi. 

Using ( 16. 4p . we can infer for X the relation 

T = ±^5^-lzi-zl 
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The discriminant A of the algebraic curve (16. ip reads in the Zi coordinates (16. 5p 

A = AzlztAlAlAl (6.6) 

with components 

^^ = zl + 2zl-25\ A2 = ^o' + 4^?-25^ 
A3 = 9zo' + ^zl{2zl - 352) + A{-2zl + 5^f . 

In particular, we have that Ai = — 2X^. 

The A-periods of the curve (16. ip with one-form (16. 3 p are taken to be 

It is not hard to explicitly evaluate the integrals in the coordinates (16. 5p for small Zj. 
The first few terms read 

9 ^i^2 . 9 dd . 9 R^2 , 

Jo\Zi) — 2 16 85^ 45^ 85^ 



9 ^2^2 9 d 9 A ji 9 dd 9 ^6 J2 
Note that the two periods are related via the identity 



(6.9) 



S, + \S, = ^(zo^ - 4zf){3zl + Azl - 45^) . (6.10) 

Inversion of (16. 9p yields the so-called mirror maps Zi{Si). 

Similarly, it is not hard to explicitly evaluate the B-periods 

no = ^/ ni= / cj, 

where A — 00 is a cutoff. Note the additional factor of 1/2 we introduced for Hq. Its 
origin can be most easily seen at hand of figure [H The period Hq (without the 1/2) in 
the quotient space corresponds in the covering space actually to 2110 and not just Hq. 
We obtain in terms of the fiat coordinates Si for the B-periods 

no(5,) =|(A2 - 25^)^^ - \ogg - 2(5i + ]^S,) log A 

+ 2(5i-^5o)log5 + Po(5o,5i), 



Ii,{S,) =|(52 _ A2)2 - 4(Si + ^S-o) log A + 2(^1 + 5o) log5 + S, log 

+ Pl(5o,5i) 



(6.11) 
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with Pi{So, Si) = Si log Si + Xln,m>o ( j2)nC-i SqSY" and Q(n, m) constants . 
Invoking the usual special geometry relation 

a5,-F(°) = n,(5o,^i)) , 

the prepotential J-'^^^ can be determined, and indeed matches the results of section 
15.21 Note that the Pi{So, Si) can be expressed in terms of the flat-coordinate Yukawa 
couplings CsiSjSk &s Pi{So, Si) = f dSidSiCsiS^Sr Closed expressions for the Yukawa 
couplings in Zi coordinates, i.e., Cz.zjz,, '■= Dz,DzjDzk^^^\zi) with Dz^ denoting the 
covariant derivative (c/., [33]), can be found to be 

^o(9^o' + 64(-5 + zl) + 8(1 - 2zl)\-l + zl) + 4^2(7 - 82? + Az\)) 



n 

'^ZQZQZO 



CziZlZl 



CzQZQZl 



n 

^ZQZIZI 



3212 

^1(34 + 8(-l + 2^2)3 + 24(-7 + 10^2) + 4^2(5 _ ig^2 ^ ^2zf) 



4X2 



Zizii-2 + 3zl){-2 + zl + Azf) 



8X2 

^o^?(4- 4(1 -2.2)2) 



4X2 



(6.12) 

The remaining couplings follow by symmetry. We have also set for simplicity g = 6 = 1. 
6.2 One-loop 

Having the tree-level data at hand, it is straight-forward to evaluate the solution to 
the 1-loop holomorphic anomaly equation of [51] 



jra)(^;/3) = ilogdetG + a^^^.;/?), (6.13) 

with Gij := ds^Zj and ; P) denoting the 1-loop holomorphic ambiguity. The 

ambiguity can be parameterized in terms of the discriminant loci Aj given in (16. 7p as 

a^^^ {z; /3) = uq log zq + ui log zi + ki log Ai + ^2 log A2 + log A3 . 

Prom the eigenvalue ensemble results of section 15.21 we deduce that under fixing pa- 
rameters Ui and Ki in the case to 



fi;i = z/o, K2 = -iyi + -, Ks = ui, 



(6.14) 
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and for Djy to 



the previous results can be reproduced. We observe that the ambiguities for both cases 
differ in general only in uq. Further, note that only in the special case of /3 = 1 and in 
the Nekrasov-Shatashvili limit (/3 — > 0) the uq of both cases are equal and the relation 

77(1) _ 

holds, as expected from the relations f l3.19p and f l5.13p . For general /3 this equality will 
not hold anymore. Also note that in the Nekrasov-Shatashvili limit (15.1 2p the 1-loop 
amplitude (I6.13P becomes purely holomorphic, i.e., 

^^'^ (^) = ^ ^0 + ^ log ^1 + ^ log (Ai A2) + ^ log A3 . 
The reason being that in this limit the 1-loop anomaly equation reduces to 

didjW^^\z,z) = 0. 

To conclude this section, it is interesting to compare the uq coefficients we found to 
the corresponding "^^g^ coefficients of section [HI given in (13.12P and ( ]3.20p . Up to an 
addition of 1/2 they match. We attribute the additional 1/2 to an artifact of our 
expansion at the 1-loop level and/or to the chosen parameterization. Mainly because 
the same mismatch by 1/2 occurs for ui, which should be equal to ^I/^''. 



7 Conclusion 

In this work we initiated the study of /3-ensembles with and Dn measure beyond 
tree- level. For that purpose, we generalized the calculation of the /3-deformed An 
partition function of [18], which makes use of a saddle-point approximation and Ward 
identities, to the Bn and Dn cases. At hand of the quartic, we found that the resulting 
free energies possess an expansion into even powers of Qs only, under a specific choice 
of 't Hooft parameters. This is as expected, since the sector is a closed period and 
should be removable via an appropriate shift of parameters following [35]. The absence 
of an odd sector in Qs allowed us to invoke the usual holomorphic anomaly equation to 
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reproduce the 1-loop sector (g^) of the quartic, albeit with new boundary conditions 
(holomorphic ambiguity) which have not appeared (to our knowledge) before. The 
boundary conditions are related to a large expansion of the Macdonald integral. 
We expect that the higher genus coefficients of the Macdonald integral expansions will 
provide boundary conditions for the higher loop amplitudes expanded near some of the 
other points in moduli space. However, so far we have not pushed the holomorphic 
anomaly calculation for the quartic beyond genus one, and it would be interesting to 
do so. Our results indicate that the /3-deformation of the Bn and Dm ensembles, which 
for /3 = 1 correspond to topological string orientifolds, is consistent. 

The Gaussian integrals also allowed us to extrapolate the and ensembles 
to toric settings. We found that the case (under an appropriate shift) should 
correspond to the topological string with boundary conditions at the conifold point 
provided by the orbifold branch of the c = 1 string. The Bn case corresponds to 
the usual circle branch, similar to the standard refined topological string related to 
An ensemble. In the toric setting, our results indicate that the equivalence between 
Bn/Dn and topological string orientifolds is not general. Rather the c = 1 moduli 
space appears to yield an independent deformation space of topological string theories, 
with only accidental correspondence to topological string orientifolds for specific simple 
geometries. 

It seems likely that one may also explore ensembles with Eq, Ej and E^ type measure 
in a similar fashion. Presumably, these ensembles are related to the three discrete points 
in c = 1 moduli space (and should not possess a consistent /3-deformation). 
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A Free energies 



The explicit results for the perturbative expansion of the free energies of the Bn and 
Dn eigenvalue ensemble with quartic potential discussed in section I5l2] are given below. 
Note that the expressions given do not take the normalization of (14. ip into account. 
The additional contribution from the normalization to the shifted free energy is given 
in terms of an asymptotic expansion of the respective Macdonald integral given in 
(E2D. 



:Fb(N,.I3) = - (8JV,?/3^ - 2Jvg/3(2/3(8JVi + 1) - 5) + ATq (8/3=(]Vi - 3)JVi + /3(8JVi - 2) + s) + (2/3 - l)JVi (2/3(]Vi - 2) + 3)) 

2 

H ^ f36iV^/3^ - 4ivg/3^(2/3(28iVi + 5) - 19) + fl2^^ (l2Nf - 16iVi + l) + ,8(24iVi - 50) + 53) 

+ JVo (-8/3^iVi (2nI - 15iVi + 19) + (^~48iVi + 200iVi + e) - ;3{88iVi + 15) + 12) - (2/3 - l)JVi (^4/3^ (/Vj - 5JVi + 5) 

+/3(15/Vi ~ 29) + 12)) 
3 

('864iV^/3* - 8iV*^^(4^(233iVi + 49) - 331) + SN^,/!'^ ( 8/3^ ( 131iV? - 102iVi + 16] - 4;9(29iVi + 103) + 383^ 

1293512^3 yuu uyyi J J 

-iN^P [ifi^ (l65ivj - 888iVi - 80) + ifi^ (ihlN^ - 558Nf + 640/Vi + I5) + 2^(965iVi + 286) - 393) + Ng (8/3^ (419ivJ 

-798iVi + 48) + 32^*" /Vi (sN^ - SS/V^ + 206/Vi - I65) + 8/3^ (112/V1 - lOSON^ + 1208/Vi - 15^ + 12/3(88iVi - 43) + 297^ 

+ (2/3 - 1)S (4/3^ (eeiVj - 314JVi + 307) + 8/9^ (siVj - 44iVi + IO8JV1 - 74^ + 6/3(86JVi - 163) + 297)) 
4 

- ^ 4^i6g4 (-6048iVo/3^ + 16JV^/3*0(4232JVi + 982) - 1549) - 4JVg3^ (4^^ (7088JVi - 3492JVi + 1033) - 23(3596JVi 
+5973) + 10325) + 8JVo/3^ (/3^ (3232iVi - 29640iVi - 4326) + 16^^ (456iVi - 1088iVi + 1275iVi + 65) + ^(18132iVi 
+6853) - 4350) - Ng ^ ^4/3^ ^27696JVi - 29580JVi + 598l) + 32/3^ (740JvJ - 7254JVi + 6IIIJV1 - 325) + 16/3* ^628-^* 

-4216iVj + 11172iVj - 7232JVi + 105) + 48/3(131JVi - 581) + 1469l) + JVq (^48/3^ (515iVj - 1808iVi - II2) 
+ 16,3^ (l474iVi - 7593iVi + 9782JVi + 195) + 32^^JVi (l4JVj' - 279iVi + 1522iVi - 2956iVi + 1769) + 8/3'' (^ASSN^ 

-7680JVj + 22104iVj - 17872JVi - IO5) + 3/3(7520JVi + 1743) - 2448) + (2/3 - l)JVi ^6/3^ ^632JVj - 2905JVi + 2799) 
+4/3^ (^279Nf - 2318JVi + 5492iVi - 3695) + 16/9'' (7JVj' - 93Nf + 398JVi - 658JVi + 353) + /3(5229JVi - 9795) + 2448)) 

2^ (93312nIi3'^ - 288JvJ;,9^(2/3(2252JVi + 559) - 1685) + 16iV^/3* ('24,9^ ( 7711JV? - 2423JVi + 125o) 

-2/3(31728JVi + 40813) + 67259) - 40iVQ/9^ ^4/9^ (l605JVi - 47267JVi - 9059) + 12,9^ (4784iVi - 8246JVi + 10442JVi 
+803) + 2/3(62029JVi + 26899) - 32463) + 4N^ 13^ ^40/3^ (35489iVi - 20156JVi + 10312) + 40^^ (5544ivJ - 66014iV^ 
+398I8JV1 - 5027) + 48^'' (3670iVi - 16580iVi + 41550JVi - 21780JVi + 863) - 30,3(4710JVi + 14933) + 223489) 
-2iVo,9 ^60/3^ ^10356JVi - 51699JVi - 6938) + 20/3^ (54900Jvj' - 186600JVi + 252646JVi + 13849) + 40/3* (4038iVi 
-61004iVj + 137454JVJ - 109848iVi - 2589) + 16^^ (2568iVj - 26610iVj + 120040JVj' - 192810iVj + 109832JVi + 945) + 
^(982802iVi + 373894) - 164865) + Nq (413^ (246049JVi - 415992JVi + 38745) + 120/3^ (4782jvf - 34959JV^ + 35464iVi 
-1034) + 8/3* (35850JV* - 334320-^]^ + 989930JvJ - 715820JVi + 7767) + 64/3''jVi (^42N^ - 1158JVj' + 9500-^]^ - 31960ivJ 
+45679JVi - 22355) + 16/3^ (^2064iVi - 47710-^* + 240920jvj' - 460650JVi + 265116JVi - 945) + 30/9(6880iVi - 4111) 
+50139) + (2^ - l)JVi (4/3^ (^27945JV5^ - 125282JVi + 119563) + 8/3^ (sSdSN^ - 46950Jvj^ + 108328iVi - 71927) 
+8/3* (ll58iV* - 14460iv]' + 59320JVi - 95329JVi + 50293) + 32^^ (21iVi - 386JV* + 2480JVj' - 7080JVi + 9025iVi 
-4081) + 6^(20555iVi - 38242) + 50139)) 
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-^d(JV,./3) =— ^ (-9iV^/i^ +iVi(3 + 2(-l + ]Vi)/3) +2]V^/3(-5 + 8(l + 2JVi)/3) - iVo{3 + 2{-S + 4iVi )/5 + 8(1 + iVi + N^)fi^)) 
2 

+ , (36Nol3^ - 4N^,(i^(-19 + 28(1 + 2]Vi)3) + JVi(12 + 15(-1 + JVi);9 + 2(3 - 5iVi + 2JvJ)/3^) 

+Ngl3{53 + 2(-77 + 12]Vi)/3 + 4(29 + 36iVi + 36iVi )/3^) - iVo(-12 + (53 + 88iVi)/3 + 2(-39 - 68iVi + 24JV^),8^ 
+8(5 + lliVi + 3iVi + 2Jvii')/3^) 
3 

H „ (-S64N^0'^ + 8JVo/3^(466/3(2JVi + 1) - 331) - SNj^ff'^ ("4,3^ (262N^ + 262JVi + ISO) - 4/3(29]Vi + 260) + 3S3) 

+4Ngl3 (20/3^ (33JV1 - 174JV1 - 109) + 4/3^ (^152Jvf + 228JVi + 622JVi + 273^ + 2/3(965JVi + 781) - 393) 

-No (419iVi + 98iVi + 398) + 8/3^ (ll2ivii' - 546N^ - 412JVi - 38l) + 32^" (bN^ + lOJvf + 80iVi + 75iVi + 37) 

+ 12/3(88iVi - 131) + 297) + JVi ^24,9^ (H-'V? - 27JVi + le) + 8/9^ (sivij' - 22jvj^ + 32JVi - I5) + 516;9(JVi - 1) + 297)) 
4 

H (eOiSNgll^ - 16iVQ ;3'*(2116/3(2iVi + 1) - 1549) + 4Ng (4/3^ (7088iVj + 7088iVi + 477l) - 2;3(3596iVi + 13449) 

+10325) - 8N^I3^ (,8^ (3232JVi - 34536iVi - 2199o) + 24,3^ (304]Vi + 456iVi + 1054JVi + 45l) + 3^(6044]Vi + 5433) 
-4350) + NgP (12,3^ (9232Jvj^ + 7404]Vi + 11455) + 32/3^ (jiON^ - 5286iVj - 5271JVi - 401o) + 16^* (628JV5' + 1256Jvf 
+ 7236iV^ + 6608iVi + 3085) + 48/3(131iVi - 1483) + 1469l) - Nq ^48^^ (515JV^ - 1296iVi - 758) + 16;3-'* (l474jvj' 
-789iVi + 5814JVi + 3023) + 8/3* (|488iVi - 3552jvf - 2496ivf - 9680iVi - 440l) + 32,8-°' (l4JVJ' + 35JVj + 490iVi 
+ 700JVJ + 937JV1 + 353) + 3/3(7520iVi + 4897) - 2448) + Ni (48/9^ (79JvJ - 191iVi + II2) + 12/3^ (93ivf - 398ivJ 
+565JVi - 260) + 8,3*' (l4JV;' - 93iVj' + 234JVi - 260iVi + lOs) + 5229^((JVi - 1) + 2448)) 

H ^° f-93312JV;^'^ + 288iV;;/3^(2252,3(2JVi + 1) - 1685) - 16iV^,3'* (s/l^ {'23133iV^ + 23133iVi + 14833) 

-2^(31728]Vi + 85145) + 67259) + 40iV^,3^ ^4/3^ (l605ivj - 60559iVi - 38558) + 4,3^ (l4352jvf + 21528ivJ + 

44618iVi + 18721) + 2/3(62029JVi + 58538) - 32463) - 4JVo;3^ ^20^^ (70978iVi + 79294JVi + 96033) + 80/3'^ (2772JVi 

-29003iVj - 31778iVi - 21998) + 32,8*' (5505iV]' + llOlOiVj + 52315JVj + 46810iVi + 20967) - 30;3(4710iVi + 34231) 

+223489) + 2JVo/3 ^60,8^ (l0356iVi - 51393iVi - 36176) + 20/3^ (54900iVi + 22896iVi + 252052JVi + 140723) 

+40/3*^ (4038JV;' - 38428ivf - 45138JvJ - 110094JVi - 50559) + 16,3'^ (2568iV^ + 6420JV;' + 63600Jvf + 88980ivJ 

+ 110042JVi + 40405) + ;3(982802JVi + 915556) - 164865) - Ng (^0^ (984196iVi + 99616iVi + 937156) + 120,3^ (4782jvf 

-19317Jvf - 10208iVi - 12766) + 8/3** (|35850JVi - 57660ivf + 373880ivf + 274580iVi + 194147) + 16/3^ (^2064JVj 

-20630iVj' - 9820ivf - 140460iVj - 115344JVi - 58455) + 64/3*' (42JVj + 126iVj + 2730iVj' + 5250iVj + 13309iVj 

+ 10705iVi + 4081) + 30/3(6880iVi - 10991) + 50139) + iVi (540,9^ (207JV1 - 494JVi + 287) + 120/3^ (393ivf - 1643iVi 

+ 2284JVi - 1034) + 24/3* (^386JVj - 2480ivf + 6050JVj - 6545JVi + 2589) + 16/3^' (|42JVj - 386iVj + 1450iVj - 2750iVj 

+ 2589iVi - 945) + 123330/3(JVi - 1) + 50139)) 

+ 0(9»). 



30 



References 



R. Dijkgraaf and C. Vafa, "Matrix models, topological strings, and supersymmetric 

gauge theories," Nucl. Phys. B 644 (2002) 3 |arXiv:hep-tli70206255 1 . 

R. Dijkgraaf and C. Vafa, "Toda Theories, Matrix Models, Topological Strings, and 



N=2 Gauge Systems," arXiv:0909.2453 [hep-th 



M. Aganagic, M. C. N. Cheng, R. Dijkgraaf, D. Krefl and C. Vafa, "Quantum Geometry 



of Refined Topological Strings," |arXiv:1105.0630 [hep-th 
L. F. Alday, D. Gaiotto and Y. Tachikawa, "Liouville Correlation Functions 
from Four-dimensional Gauge Theories," Lett. Math. Phys. 91 (2010) 167 



arXiv:0906.3219 [hep-th 



N. A. Nekrasov and S. L. Shatashvili, "Quantization of Integrable Systems and Four 
Dimensional Gauge Theories," larXiv:0908.4052 [hep-th] 



M. Aganagic and S. Shakirov, "Knot Homology from Refined Chern-Simons Theory," 
arXiv: 1105.51 17 [hep-th] [ 

D. Ghoshal and C. Vafa, "C = 1 string as the topological theory of the conifold," Nucl. 
Phys. B 453, 121 (1995) [arXiv:hep-th/9506122| . 

S. K. Ashok, R. Corrado, N. Halmagyi, K. D. Kennaway and C. Romelsberger, "Unori- 
ented strings, loop equations, and N = 1 superpotentials from matrix models," Phys. 
Rev. D 67 (2003) 086004 |arXiv:hep-th/0211291| . 

K. A. Intriligator, P. Kraus, A. V. Ryzhov, M. Shigemori and C. Vafa, "On low rank 
classical groups in string theory, gauge theory and matrix models," Nucl. Phys. B 682 



(2004) 45 arXiv:hep-th/0311181 



N. Nekrasov and S. Shadchin, "ABCD of instantons," Commun. Math. Phys. 252 (2004) 
359 | arXiv:hep-th70404225] . 

M. Aganagic and K. Schaeffer, "Orientifolds and the Refined Topological String," 
[arXiv: 1202. 4456 [hep-th 



K. Landsteiner, C. I. Lazaroiu and R. Tatar, "(Anti)symmetric matter and superpoten- 
tials from IIB orientifolds," JHEP 0311 (2003) 044 |arXiv:hep-th/0306236] . 
K. Landsteiner and C. I. Lazaroiu, "On Sp(0) factors and orientifolds," Phys. Lett. B 



588 (2004) 210 |arXiv:hep-th/0310111] . 

A. Brini, M. Marino and S. Stevan, "The Uses of the refined matrix model recursion," 



arXiv:1010.1210 [hep-th] 



V. Bouchard, A. Klemm, M. Marino and S. Pasquetti, "Remodeling the B-model," 
Commun. Math. Phys. 287 (2009) 117 [arXiv: 0709. 1453 [hep -th^. 

M. Marino, "Open string amplitudes and large order behavior in topological string 
theory," JHEP 0803 (2008) 060 |arXiv:hep-th/0612127] . 



31 



[17] A. S. Alexandrov, A. Mironov and A. Morozov, "Partition functions of matrix mod- 
els as the first special functions of string theory. 1. Finite size Hermitean one ma- 
trix model," Int. J. Mod. Phys. A 19 (2004) 4127 [Teor. Mat. Fiz. 142 (2005) 419] 
[arXiv:hep-th /0310113l. 

A. Morozov and S. Shakirov, "The matrix model version of AGT conjecture and CIV-DV 



prepotential," JHEP 1008 (2010) 066 arXiv:1004.2917 [hep-th 



M. Aganagic, A. Klemm, M. Marino and C. Vafa, "Matrix model as a mirror of Chern- 
Simons theory," JHEP 0402 (2004) 010 | arXiv:hep-th /0211098l. 

A. Klemm, M. Marino and S. Theisen, "Gravitational corrections in supersymmetric 
gauge theory and matrix models," JHEP 0303 (2003) 051 jarXiv:hep-th/0211216] . 
I. G. Macdonald, "Some Conjectures for Root Systems," SIAM J. Math. Anal., Vol. 13, 
No. 6, November 1982 

E. M. Opdam, "Some applications of hyper geometric shift operators," Invent. Math. 98, 
275-282, 1989 

E. M. Opdam, "Dunkl operators, Bessel functions and the discriminant of a finite Co- 
exter group," Compositio Math. 85, 333-373, 1993 

M. Aganagic, R. Dijkgraaf, A. Klemm, M. Marino and C. Vafa, "Topological strings and 
integrable hierarchies," Commun. Math. Phys. 261 (2006) 451 |arXiv:hep-th/0312085| . 
J. Harer and D. Zagier, "The Euler Characteristic of the Moduli Space of Curves," Inv. 
Math. 85 (1986) 457-485 

L. Chekhov and A. Zabrodin, "A Critical matrix model for nonoriented string," Mod. 
Phys. Lett. A 6 (1991) 3143. 

I. P. Goulden, J. L. Harer and D. M. Jackson, "A geometric parametrization for the 
virtual Euler characteristic of the moduli space of real and complex algebraic curves," 
Trans. Am. Math. Soc. 353, 4405 (2001) 

H. Ooguri and C. Vafa, "World sheet derivation of a large N duality," Nucl. Phys. B 



641 (2002) 3 [arXiv:hep-th702052 97 1 . 

D. Krefl, S. Pasquetti and J. Walcher, "The Real Topological Vertex at Work," Nucl. 



Phys. B 833 (2010) 153 [arXiv:0909.1324 [hep-th 



D. Krefl, "Wall Crossing Phenomenology of Orientifolds," larXiv: 1001. 5031 [hep-th] 



D. J. Gross and I. R. Klebanov, "One-dimensional String Theory On A Circle," Nucl. 
Phys. B 344 (1990) 475. 

M. -X. Huang, A. Klemm and S. Quackenbush, "Topological string theory on compact 
Calabi-Yau: Modularity and boundary conditions," Lect. Notes Phys. 757 (2009) 45 
[arXiv:hep-t h/0612125| . 

M. Bershadsky, S. Cecotti, H. Ooguri and C. Vafa, "Kodaira-Spencer theory of gravity 
and exact results for quantum string amplitudes," Commun. Math. Phys. 165 (1994) 
311 [arXiv:hep-th/9309140] . 



32 



D. Krefl and J. Walcher, "Extended Holomorphic Anomaly in Gauge Theory," Lett. 
Math. Phys. 95 (2011) 67 |arXiv: 1007.0263 [hep-th]] . 

D. Krefl and J. Walcher, "Shift versus Extension in Refined Partition Functions," 



arXiv:1010.2635 [hep-th 



M. -X. Huang and A. Klemm, "Direct integration for general O backgrounds," 



arXiv: 1009. 1126 [hep-th 



J. Walcher, "Extended holomorphic anomaly and loop amplitudes in open topological 
string," Nucl. Phys. B 817 (2009) 167 |arXiv:0705.409 8 [hep-th]|. 

D. Krefl and J. Walcher, "The Real Topological String on a local Calabi-Yau," 



arXiv:0902.0616 [hep-th 



P. H. Ginsparg, "Curiosities at c = 1," Nucl. Phys. B 295 (1988) 153. 

R. Gopakumar and C. Vafa, "Topological gravity as large N topological gauge theory," 

Adv. Theor. Math. Phys. 2 (1998) 413 |arXiv:hep-th/9802016 . 



D. Ghoshal, D. P. Jatkar and S. Mukhi, "Kleinian singularities and the ground ring of 

C=l string theory," Nucl. Phys. B 395 (1993) 144 |arXi v:hep-th/9206080j . 

D. Krefl and S. -Y. D. Shih, "Holomorphic Anomaly in Gauge Theory on ALE space," 



arXiv:1112.2718 [hep-th 



B. S. Acharya, M. Aganagic, K. Hori and C. Vafa, "Orientifolds, mirror symmetry and 



superpotentials," |arXiv:hep-th/0202208 

J. Walcher, "Evidence for Tadpole Cancellation in the Topological String," 



|arXiv:0712.2775 [hep-th 



A. Mironov, A. Morozov, A. Popolitov and S. .Shakirov, "Resolvents and Seiberg-Witten 



representation for Gaussian beta-ensemble," |arXiv:1103.5470 [hep-th 



I. P. Goulden and D. M. Jackson "Maps in Locally Orientable Surfaces and Integrals 
over Real Symmetric Surfaces," Can. J. Math. Vol. 49(5), 1997 pp. 865-882 
H. Ita, H. Nieder and Y. Oz, "Perturbative computation of glueball superpotentials for 
SO(N) and USp(N)," JHEP 0301 (2003) 018 [hep-th/0211261]. 

R. A. Janik and N. A. Obers, "SO(N) superpotential, Seiberg-Witten curves and loop 
equations," Phys. Lett. B 553 (2003) 309 [hep-th/0212069]. 

H. Fuji and Y. Ookouchi, "Confining phase superpotentials for SO / Sp gauge theories 
via geometric transition," JHEP 0302 (2003) 028 |arXiv:hep-tIi70205301] . 
R. Dijkgraaf, S. Gukov, V. A. Kazakov and C. Vafa, "Perturbative analysis of gauged 
matrix models," Phys. Rev. D 68 (2003) 045007 |arXiv:hep-th/0210238| . 
M. Bershadsky, S. Cecotti, H. Ooguri and C. Vafa, "Holomorphic anomalies in topolog- 
ical field theories," Nucl. Phys. B 405 (1993) 279 |arXiv:hep-th/9302103| . 



33 



